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Abstract 

We present first-principles calculations on bulk CdSe and CdSe nanowires with diameters of up to 22 A. 
Density functional linear combination of atomic orbitals and plane wave calculations of the electronic and 
structural properties are presented and discussed. We use an iterative, symmetry-based method to relax the 
structures into the ground state. We find that the band gap depends on surface termination. Vibrational prop- 
erties in the whole Brillouin zone of bulk CdSe and the zone-center vibrations of nanowires are calculated 
and analyzed. We find strongly size-dependent and nearly constant modes, depending on the displacement 
directions. A comparison with available experimental Raman data is be given. 

PACS numbers: 61.72.uj, 63.22.-m, 63.22.Gh 
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I. INTRODUCTION 



Nanocrystals and nanowires (NWs) have become an important field in solid state physics, es- 
pecially due to the recent advances in the growth methods. For example, the fabrication of CdSe 
NWs or nanorods with small diameter deviations and diameters down to a few nanometers has 
been accomplished.-^ On a nanometer scale confinement effects strongly influence the electronic 
structure. While there has been considerable effort in understanding the electronic and optical 
properties from theory^^ and experiments^ '^ there have been little efforts to understand the vi- 
brational properties of such nano structures and of bulk CdSe, either. The phonon dispersion of 
wurtzite CdSe is only partly available from experiments^. To our knowledge there is no first- 
principles calculation of the phonon dispersion available for hexagonal CdSe. Only a bond-charge 
model calculation with nearest neighbor interactions is available.- Recent studies on hexagonal 
materials, like, e.g. graphite, revealed that for a proper description of the vibrational properties of 
hexagonal systems interactions up to four nearest neighbors are needed.- Confinement effects have 
been discussed for nanospheres.-. Recently, Raman data on small axial-symmetric CdSe nanorods 
have been presented.— 

In the present work we concentrate on the vibrational properties of CdSe bulk and NWs. We 
compare nanowires with diameters between 5 and 22 A. Structural and electronic properties de- 
pending on the diameter and surface termination are investigated. We show how the counterparts 
of the optical bulk modes evolve in frequency with decreasing diameter. We find phonon modes 
with higher frequencies than in the bulk, that vibrate mainly on the surface. 

II. CALCULATIONAL DETAILS 

Two different DFT codes were used for the calculation of bulk properties and NW properties: 
for the phonon dispersion of bulk material we used linear response, which is implemented in 
ABINIT.— i^^ii^ For large unit cells, as required by one-dimensional systems, however, plane 
wave basis sets are computationally very expensive. Therefore we used for NWs the SIESTA 
code that uses atom-centered confined numerical basis functions .-^^^ 

The NWs were generated by cutting appropriate discs out of a wurtzite crystal. The cross 
section of the relaxed unit cells is shown in Fig. \T\ We used the local density approximation^^. 
Using gradient approximation should not influence the tendencies will impose. Pseudopotentials 
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FIG. 1: (Color online) Unit cells of all calculated hydrogen-coated NWs. The number of CdSe-pairs per 
unit cell is indicated. As can be seen, near the surface, the Cd atoms (magenta, dai^k) move more inwards, 
as compared to bulk region.— 

were generated with the TrouUier-Martins scheme^ for the following valence-state configurations: 
Se 45^(1.89) 4p'^(1.89) Cd 55^(2.18) 4(i^°(2.5), where the value in parenthesis indicates the 
pseudopotential core radii in bohr. The valence electrons were described by a double-^ basis 
set plus an additional polarizing orbital. The localization of the basis followed the standard split 
scheme and was controlled by an internal SIESTA parameter, the energy shift, for which a value 
of 50meV was used. This resulted in basis functions with a maximal extension of 3.45 A(Se), 
4. 13 A(Cd) and 3.2 A(H). Periodic images of the NWs were separated by at least 20 A. Real space 
integrations were performed on a grid with a fineness of 0.2 A, which can represent plane waves 
up to an energy of 80 Ry. For bulk CdSe a (6 x 6 x 4) Monkhorst-Pack-° mesh in reciprocal space 
was used, whereas for wires a minimum of 16 /c-points equally spaced along the ID Brillouin zone 
was used. The phonon calculations were performed with the method of finite-differences.^' For 
the calculation of the LO splitting we used the experimental dielectric high-frequency constant-^ 
etL = 6.0 and a modified version of SIESTA described in Ref 23 



The lattice parameters for wurtzite CdSe obtained with ABINIT were a=4. 29 A and c=7. 00 A; 
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SIESTA yields a=4.34A and c=7.09A, both in good agreement with the experimental values 
a=4.2999 A and c=7.0109 A from inelastic neutron scattering (INS).- We also calculated the bulk 
modulus and obtained the values 55.58 GPa (ABINIT) and 55.54 GPa (SIESTA) again in good 
agreement with the experimental value 53.4 GPa.—. 

Symmetry considerations helped to reduce the calculation time. NWs can be described in terms 
of line groups,— which are the one-dimensional analogon to space groups. The NWs from Fig. [U 
with a hexagon in the center (upper row) belong to the line group L(6)^mc which has the generator 
{C6|t/2}. This manifests itself in a screw axis. The NWs from Fig. [T] with three 120° bonds in the 
center (lower row) belong to the line group L3m which has the generator {C3IO}. In large systems 
internal stresses can occur, therefore we symmetry based, iterative the relaxation process: After 
the NW was relaxed, we generated the whole wire by applying symmetry operations to a minimal 
set of atoms in the unit cell. This symmetry generated cell was relaxed again. This iterative 
process was repeated until all forces were below 0.04 eV/A after generating the wire. Although 
the relaxation process is considered to not change the symmetry we observed better results with 
the generating method: Frequencies of phonons belonging to two-dimensional representations 
(degenerate frequencies) were separated by less than 0.05 cm"^. 



III. RESULTS AND DISCUSSION 



The unit cell length c of the NWs changes with diameter and with the addition of a hydrogen 
passivation shell around them. We plotted the unit cell length c over the inverse diameter 1/d in 
Fig.[2l The diameter is defined by the average distance of the outermost Se atoms to the rotational 
axis. The bare wires have a larger unit cell parameter c than the passivated NWs. The largest 
increase of 3.5 % was found for the (CdSe)i3 wire. For the largest wire, the (CdSe)54, the increase 
is 1.1%. We find a linear relationship, which is best described by the parameters Cnw('^) = c — m/d 
with the bulk lattice constant c = 7.09 A and m = 3.25 A^ (passivated) and m = 1.74 A^ (bare) 
NWs. 

The most dramatic changes in the atomic structure happen near the surface. The surface of all 
calculated NWs undergo a reconstruction. This reconstruction was first reported by Ref. [l7| via 
extended x-ray absorption fine structure spectroscopy and first-principles calculations. During this 
reconstruction the Cd-atoms rotate into the surface to lower the energy as can be seen in Fig. [U 
The passivation with hydrogen atoms does not qualitatively change this reconstruction. After 
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relaxation, the Cd-Se bonds show a broad distribution in lengths. In the core the bond lengths 
deviate from the bulk value by < 1 % , while those on the surface change by up to 5 % . The 
change in bond length causes changes in the force constants which will be discussed later. 

CdSe has a direct electronic band gap. It is known from absorption or photoluminescence mea- 
surements that the band gap of CdSe nanocrystallites and NWs shows a strong size dependence.- 
In Fig. |2] (lower) we plot the band gap increase (compared to the calculated bulk value) over the 
inverse diameter. As can be seen the passivation influences the electronic band structure; it gener- 
ally increases the band gap at the T point. At first sight this is a strange behavior as the H-atoms 
add electronic states to the unit cell. However, at the same time the lattice parameter c decreases 
with the addition of the H-atoms which is a stronger effect. Thus the lateral confinement increases 
and leads to a net increase of the band gap. The data from the wires with passivation is in excellent 
agreement with the shown semiempirical pseudopotential calculation (SPC) data from Ref. [ij that 
are based on experimental data. 

The phonon dispersion calculated with ABINIT is shown in Fig. |3l Also shown are ex- 
perimental frequencies obtained with INS along the F— A direction which was multiplied by 



a/ (116/112.4) to account for the use of ^^^Cd isotope instead of natural Cd with an atomic weight 
of 112.4U.- The F-point phonons calculated with SIESTA are also plotted; we compare them 
later with the NW frequencies. All frequencies were multiplied by a constant factor to match the 
experimental bulk LO frequency at the F-point. The dispersion along F — K and F — M shows 
an overbending near the F-point, which can be found in other wurtzite crystals as well.-^^ When 
treating the vibrations in the zone-folding approximations, it is important to know whether or not 
there is overbending.— 

We have plotted a normalized histogram of the F-phonons of a representative NW into the bulk 
phonon density of states (DOS) in Fig. [3l Acoustic modes below 75cm^^are shifted to higher 
frequencies, whereas the optical modes around 200cm~^are red-shifted, a general trend, e.g., as 
reported by Sun in Ref 129. Additionally, optical modes in NWs are found with frequencies higher 
than the LO-mode in the bulk: These modes vibrate mainly on the surface, their frequency exceed- 
ing the bulk LO by at most 15 cm^^ The reason for the increase lies in the higher coordination of 
the reconstructed Cd atoms as an inspection of the force constants for the outermost Cd atoms for 
the slab shows. The sum of the forces on this atom when slightly displaced is higher by 4 % than 
that of an Cd atom in in the center for equal displacement. Force constants of atoms lying in the 
2nd layer or closer to the rotational axis do not significantly differ from those in the bulk. This is 
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similar to the subsurface modes predicted for the zirconium surface.— 

We now focus on optical phonon modes that exist in bulk material and their counterparts in 
NWs. The evolution in phonon frequency with changing NW diameter is shown in Fig.lH where 
we plot the phonon-frequency difference between the NW and bulk over the inverse diameter. All 
considered bulk modes that belong to one-dimensional representations show a strong systematic 
redshift for reduced diameters (upper plot).— In contrast, the modes belonging to two-dimensional 
representations show a smaller size dependence (lower plot). The rigid shift of the i^i-mode 
will be explained later. This different size dependence has also been predicted for Si NWs in 
Ref. [sil They find the LO mode with eigenvectors parallel to the NW axis shows a stronger size 
dependence, than the two TO modes with eigenvectors perpendicular to the axis. This is the same 
here, all modes belonging to one-dimensional representations have eigenvectors almost parallel to 
the wire axis and depend strongly on size. 

Polar semiconductors are known to have different frequencies for LO and TO modes at the 
F-point. The reason for this is that long-wavelength optical longitudinal atomic vibration induce 
a polarization field that contributes to the restoring forces due to the long-range Coulomb inter- 
action. It has been shown, that this Coulomb interaction can be treated separately from the short 
ranged dynamical calculation. In detail the dynamical matrix can be split into an analytical and 
a non-analytical part, the latter containing the long-range Coulomb forces. Its evaluation is per- 
formed as an Ewald summation. ■^^'^^ This non-analytical part depends on the unit cell volume and 
the dielectric high-frequency constant eoo- In bulk CdSe with an experimental dielectric constant^ 
of eUo = 6.0 a splitting of 32cm~^ is obtained in good agreement with the experimental value of 
39 cm~^.- We apply the same model to NWs and take the wire cross section times unit cell height 
as cell volume. To account for the one-dimensionality we only sum over k- vectors along the wire 
axis. The same dielectric constant e^L is used.'^ We obtain a splitting of only a few cm~^ for larger 
wires; much smaller than in the bulk. We attribute the underestimated splitting to the fact that the 
displacements within the LO mode in the NW are not well ordered. They are parallel only in the 
NW-core and even antiparallel near the surface. Consequently the polarization as sum of dipoles 
becomes weaker compared to bulk. In Fig. |4] the energy difference of the Ai-modes is calculated 
without any contribution from the long-ranged Coulomb interactions, neither for bulk nor NWs. 
The bulk frequency reference becomes then 177 cm^^(calc. with SIESTA). For this mode exper- 
imental Raman data is available for nanorods with aspect ratios of > 5 and diameters between 4 
and 8 nm.-^^ A linear fit cUnw — ujbuik = m/d leads to m = —149 cm^^A. The calculated frequen- 



6 



cies of the Ai-modes shown in Fig. |4] are in good agreement with the extrapolated experimental 
fit. 

A rigid shift of the £'i-mode was predicted by Fuchs and Kliewer— in a calculation of an ionic 
crystal slab. They showed that phonons with a polarization perpendicular to the slab surface have 
frequencies similar to the LO-mode in the bulk. This is attributed to surface charges. A character- 
istics of this mode is that the Cd and Se sub-lattices move into opposite directions. To verify this 
we calculated the phonon frequencies of a thin slab extending into the y and ^-directions. The non- 
analyticity is taken along z. We find that the degeneracy of the Ei-mode is lifted into two modes: 
One has a polarization perpendicular to the surface and a frequency of the LO in the bulk; the other 
one has a polarization vector parallel to the surface and a frequency of the Ei-mode in bulk. The 
rigid shift of the Ei mode in NWs is thus caused by surface charges in a one-dimensional system 
and the mode is degenerated again, but now has the frequency of the LO mode in the bulk. The 
.El-modes in Fig. |4]b) show a size dependence. This can be qualitatively explained considering 
the electric field generated by the surface charges, that becomes weaker for larger diameters. 

IV. SUMMARY 

We presented first-principles calculations of CdSe bulk and NWs. We used an iterative relax- 
ation process that reproduces well the degenerate phonon modes. The effect of a H-passivation on 
NWs resulted in smaller unit cells and larger band gaps. We analyzed the vibrational properties of 
NWs. A reconstruction in the surface leads to higher force constants, resulting in surface modes 
with frequencies up to 15 cm^ higher than in the bulk. The size dependence of phonon modes that 
correspond to bulk optical F modes was analyzed. We find that modes which vibrate mainly along 
the wire axis showed a stronger size-dependence than modes with eigenvectors perpendicular to 
the axis. 
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FIG. 2: (Color online) a) Evolution of the unit cell lengths of the NWs with and without H-atoms. Bulk 
CdSe is considered as an inverse diameter of 0. The number of CdSe-pairs is indicated, b) Band gap 
difference AE = i^^w — ^buik of the NWs with and without H-passivation. Crosses denote data from 
semiempirical pseudopotential calculation (SPC) from Ref . Q. 
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FIG. 3: (Color online) Phonon dispersion curves and DOS for bulk CdSe as calculated with ABINIT 
(solid lines). Diamonds are experimental data from Ref. |5|- Circles are T-point frequencies calculated 
with SIESTA. In the bulk phonon DOS (hne) we plotted the T-DOS of the ideal (CdSe)54-wire (histogram). 
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FIG. 4: (Color online) Phonon energy difference cjnw — <j-'buik over inverse diameter for all optical bulk 
modes. On the left are the modes of a slab with a thickness of 20 A which corresponds to the diameter of 
the biggest calculated NW. a) Phonons belonging tb21 -dimensional representations, b) Phonons belonging 
to 2-dimensional representations. 



